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POTENTIAL THEORY OF RADIATION 

BY 

Huei-Huang Chiu 

Professor of Mechanical Engineering ^ 
University of Illinois at Chicago 
Chicago, Illinois 


ABSTRACT 

This study aims to develop a new theoretical method by which the structure of a 
radiation 6eld can be predicted by a radiation potential theory, similar to a classical 
potential theory. The introduction of ai scalar potential is justified on the grounds that 
the spectral intensity vector is irrotational. The vector is also solenoidal in the limits of 
a radiation field in complete radiative equilibrium or in a vacuum. This method provides 
an exact, elliptic type equation that will upgrade the accuracy and the efficiency of the 
current CFD programs required for the prediction of radiation and flow fields. 

A number of interesting results emerge from the present study. First, a steady state 
radiation field exhibits an “optically modulated inverse square law” distribution character. 
Secondly, the unsteady radiation field is structured with two conjugate scalar potentials. 
Each is governed by a Klein- Gordon equation with a frictional force and a restoring force. 
This steady potential field structure and the propagation of radiation potentials are con- 
sistent with the well known results of classical electromagnetic theory. The study also 

recommends the extension of the radiation potential theory for spray combustion and 
hypersonic flow. 
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1. INTRODUCTION 


Thermal radiation plays an important role in the broad area of engineering and sci- 
entific applications involving high temperature flow processes. In high performance rocket 
engines industrial boilers and furnaces, as well as in many astrophysical flow fields^, the 
radiative tia.isfer frequently makes a dominant contribution to the energy redistribution, 
mass, and the momentum transfer in the participating media. 

Recently, an assessment of the impact of radiation in rocket engine performance was 
conducted by this writer^ It was concluded that the radiative heat loss and the en- 
hancement of the burning rate of droplets by radiation can significantly affect the overall 
performance of low-intermediate enthalpy engines. The results of this preliminary study 
were brought to the attention of the research community of liquid propulsion technology. 
Subsequently, the JANNAF workshop^ on “Radiation effect on flow characteristic in com- 
bustion chambers” was conducted to indentify the problem areas that could have potential 
impacts on rocket engine hardware design, performance, life cycles thermal material fa- 

tigue, and reliability as well as the radiative modulation on spray flow field behavior and 
combustion instability. 

The present study is the continuation of the radiation research which aims to promote 

a state-of-the-art understanding of the phenomena and to upgrade the predictive capability 

as well as the accuracy of Computational Fluid Dynamics (CFD) codes with radiative heat 
transfer. 

The specific goals of this study are (l) to establish a theoretical link between two 
principal theories of radiation, electromagnetic (EM) theory vs. radiative transfer (RT) 
theory, and (2) to develop a unified theoretical methodology that can be incorporated in 
a computational aiogrithm and a grid structure similar to those of the participating fluid. 
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The 6rst problem treated is a fundamental issue of establishing a unified view of 
radiation by examining the structures of the radiation field predicted by two rival theories; 
the EM and RT models. 

The most remarkable findings that emerged from this study are (l) the radiation field 
intensity vector is irrotational and therefore the field' vector is derivable from a scalar 
potential, (2) a steady state field structure predicted by RT theory obeys an elliptic type 
equation in much the same manner as an electro- magnetic static field and (3) a non-steady 
radiative field is governed by a wave equation which has the same general features of an 
electro-magnetic field characterized by the Maxwell equation. 

The second problem in this study is the establishment of an exact radiation poten- 
tial equation for a radiation modulated Computational Fluid Dynamics program. This 
study is further motivated by the fact that the radiation equations currently adopted in 
the majority of the Computational Fluid Dynamics (CFD) approaches are approximate 
equations®; e.g., Two-Flux model, Six-Flux model, discrete-ordinate or zonal method. The 
accuracy of these equations in general is not compatible with that of the conservation laws 
and associated submodels which are more sophisticated in the process description and 

mathematical characterization; e.g., A - c model, combustion-turbulence interaction and 
spray-droplet models®. 

Since the equation governing the radiation potential in the present theory is of an el- 
liptic type, a numerical algorithm and grid net work similar to those of the Navier-Stokes 
equations, or turbulent flow conservation laws can be adopted for the simultaneous predic- 
tion of the radiation and flow field. These advantages facilitate computational efficiency 

and accuracy. 

The potential theory serves to provide a new physical perception and an useful 
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methodology toward the understanding of radiation processes in high temperature flow 
phenomena. It also provides a gauge for the comparative sissessment of the validity and 
limitation of the classical electromagnetic theory and the phenomenological radiative trans- 
fer approach. 

This writer recommends that the steady state potential theory be implemented for 
the prediction of the combustion flow field in a selected liquid rocket combustion chamber, 
such as TRW’s variable thrust engine, to test the viability of the present method. It is also 
suggested that this theory be extended to treat the interaction of radiation with droplets 
and particles in a spray combustion environment. 

Additionally, spontaneous or non-steady radiation processes in combustion environ- 
ments may be examined to assess the radiation induced ignition, vaporization and com- 
bustion in advanced propulsion systems. 


2. BASIC PROPERTIES OF SPECTRAL INTENSITY VECTOR AND RA- 
DIATION POTENTIAL 


The radiative transfer equation^ (RTE) which governs the spectral intensity is given 
by 

^ "b ^ /(/) (1) 

where c is the speed of light in a vacuum, lu is the spectral intensity of radiation in 
the frequency interval u and u + du, 0 is the volumetric extinction, which is the sum of 


the absorption coefficient amd scattering coefficient <Ti,, J is the unit propagation vector 
and S„ is the sum of emitted and in-scattered radiation given by 

Su ={t1u/0u) + {<r^/4TT0„)x 



v)r^>{s')dn'dv' 
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in which T]^, is the emission coefficient v !«» e 

nt, Xi/ IS the scattering function and Q 


angle. 


is a solid 


2.1 Irrotationality of Spectral Intensity Vector at Steady State 

The spectral intensity vector is defined as the vector whirK K .u 

« (.ne vector which has the magnitude equal to 

the spectral intensity and the direction coincides with that nf tK 

inciaes with that of the propagation of a ray 

of light. The spectral intensity vector at a steady state ran k u 

any state can be shown to be “irrotational” 

by the following mathematical procedure. 

By a simple quadrature the solution of Eq.(l). is given by 

Ioexp(^ 

Note that Eqn (3) is an intergral equation if 5. contains the scattering term which 
. represented by the third term on the right hand side of Eq.(2). In order to show the 
rotationality of the spectral intensity vector, the curl of I^S must vanish. 

V X (I^s) = exp( • dr'j 

- 0us 1^ exp 

+ xs = _ 0^s,)l XI = 0 

The irrotationality of a spectral intensity vector allows one to introduce a ^radiation 
potential from which the spectral intensity vector can be calculated ae follows 


(4) 


= Iu9 


Additionally, a ray of light propagates in the direction 
surface, as shown in the following. 


(5) 


normal to an equipotential 
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Since 5 IS a unit vector, the magnitude of the spectral density vector is given by 


Thus from Eq.(5), the unit propagation vector s is 

3 = ^ 

~ ( 7 ) 

Hence the propagation vector coincides with the unit vector, n, normal to an equipo- 

tential surface. 

By substituting (5) into (l) one finds that the radiation potential obeys the fol- 
lowing “Radiation Potential Equation” (RPE) 

The structure of a radiation field and spectral intensity vector can be predicted by 
the solution of Eq.(8). A general solution of Eq.{8) will be discussed in section 3. 

3-2 Solenoidal Characteristics of Spectral Intensity Vector for Two Limiting 
Cases 

An observation of E,.(l) auggeeta that when the right hand aid. of the equation 
vaniahea, the spectral intensity vector i, solenoidal or divergent free, providing that the 
radiation i, independent of time. This special case corresponds to (l)a complete radiative 
equilibrium, and (2) non-participating media i.e. 0 ^ = 0 including a vacuum space. 

In the special cases listed above, the radiation potential 4 : satisfies the Laplace equa- 

feion 

V V,. = 0 
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Thu, the determination of the dietribution of the epeetrai intensity vector for these 
limited case, i, equivalent to the boundary value problem of a “Laplace" potential field. 
Principle tool, and the result, of classical potential theories (such a, the method of image, 
conformal mapping, a, well a. Green functions for Dirichlet and Neumann boundary value 
problems) can be used with or without modification for the prediction of the fiow coupled 
or decoupled radiation 6eld problems. 

3 . POTENTIAL THEORY OF STEADY STATE RADIATION 

The determination of a radiation intensity in a domain bounded internally or exter- 
nally by a physical surface is reduced to an appropriate boundary value problem of RPE, 
Eq.(8). 

A general solution of Eq.(8), “optically modulated inverse square law", is discussed in 
the following section. 


3.1 Optically Modulated Inverse Square Law 

In order to determine the radiation field one has to predict I, simultaneously with the 
conservation laws of a flow field because the source term S„ appearing on the right hand 
side of Eq.(8) depends , in general, on the properties of a gas flow fieid and radiation. 

In the following anaiysis, however, the S„ will be treated as a source term so that the 

formal solution of Eq.(8) can be expressed in an internal form which contains 5„ 

in ter grand. 


> 1 , in an 


The general solution is first expressed by the method of potential splitting as follows 

'7d„ = exp|-^ He3'.d7'jvg.„ (g) 
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Subsequently by substituting (9) into (8) one obtains 


= exp 



df' 




(10) 


Thus a particular solution of Eq.(lO) is what may be termed as an “optically retarded 
inverse square law potential” given by 


0*. 


-m 


0uS^{r') exp(/J • dr ) 


|r-r'| 


dV 


( 11 ) 


The: opectral intensity vector is obtained by substituting Eq.(ll) into (9) as follows 

= V0, -- III dr^dV (12) 

Observation of Eq.(2) suggests that the intensity decay pattern has two primary char- 
acteristics; inverse square and exponential decay in source strength due to an opacity 
factor. This decay pattern is one of the basic characteristics of the radiation heat transfer 
and will be referred to as “optically modulated inverse square law”. 


3.2 Steady State Radiation Potential Solution 


A general solution satisfying a prescribed potential and its normal derivative at a 
physical boundary shown in Fig.l can be constructed from the solution of an associated 
potential 0^,. Since satisfies the Poisson equation, the general solution is given by the 
following classical expression 




/// 


0vSi,{r') exp(fS^ 0^1” ■ dr”) 


\f-r\ 


dV' 


-t- 



0^{r o) 1 

dn !r — fo| 


(13) 
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a3 


where 0i/{fo) and is the value of the potential and its normal derivative at fg, 

located on the boundary. The numerical values of ^ are related with those of ^ 

5 Tl 

follows 

— (r ^) = _ijm n ■ = _lim n • exp J ■ dr'j 

Thus the numerical value of ^ at r = rg is equal to that of i.e. 

drpujro) ^ d<t,^r,) 
dn dn 

The numerical values of rpi,{fo) must be calculated from the prescribed value of 
by integrating Eq.(9) as follows 


’/'u(r) — 0t,(ro_o) + f exp{ f ■ df')dti;^ 
•'^0,0 J 

where tq q is a reference point on the wall. 


(14) 


An additional step is required to determine the value of 0„(fo) in terms of <^^(fo). 
firstly, Eq.(14) will be specialized to the boundary by replacing f by fg, i.e. 


01/ (fo) = 0i/(fo,o) + f exp(A)d<^i, 


(15) 


where 


A=/ 0^1' dr‘^ (16) 

Since the distribution of a potential is prescribed as function <f>u i.e. 


01/ — 0I/(A) 

Finally, by subsituting (17) into (16) one obtains 


0i/(ro) = 0i/(ro,o) + [ e^/(A)0A 


(18) 


8 



where 


fW = 




dX r=r„ 

The spectral intensity vector is derived from Eqs.(9) and (13) 2 is follows 

1^3 = = - J j J df”^l'dV' 


(19) 


where 


s; = 


= 


X, - I. 


- x'jV 


X, — Xoi 


- x'y 


( 20 ) 

( 21 ) 


and 01, (ro) is given by Eq.(l8). 


In general, the method of Green functions can be used to predict the solution cor- 
responding to the boundary value problem with a prescribed normal derivative of the 
radiation potential. For instamce, the solution with a prescribed inhomogeneous Neumann 
boundary condition is given by 

= III VG(r/rO/?.,S^(rOexp(^- 
+ j j (r/Fo) exp^- 0^1" ■ dr"^dY^^ 


( 22 ) 


where 


<50u(»'o) 


dn 


dn 


- Ill VG{rJr')0,S^{r)dV 


(23) 
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The Green function G satisfies the following equation 


V^G(r/r') = -47T ^(r - ?') (24) 

and homogeneous Dirichlet condition i.e. G(fo/f^) = 0 

3.3 Structural Equivalence Between the Radiation Potential Theory and the 
Electro-magneto Static Theory 

The fact that the radiative potential equation Eq.(8a) and (8b) is an elliptic type 
suggests the structural similarity between the radiation field predicted by the Radiation 
Potential Theory, and the electro-magnetic static field which is governed by the Laplace 
or the Poisson equation in electromagnetic theory. However, it must be pointed out that 
the radiation potential <^u and spectral intensity vector /„5 are not the same physical 
quantities as the electro-magnetic static potential and electro-magnetic intensity. Hence, 
complete equivalence between the two theories will require the knowledge of the functional 
relationship between those of the electromagnetic properties and those of the radiation 
field. This issue will be discussed further in section 4. 

4. POTENTIAL THEORY OF NON-STEADY RADIATION 

Unsteady radiation phenomena involve the propagation and interaction of light with 
matter in a participating fluid. The objective of this section is to develop a hyperbolic type 
equation that governs what may be termed a “conjugate spectral intensity” . The proposed 
formulation of a “Radiation Wave Equation” exhibits another interesting insight into the 
structural equivalence between the radiative transfer and electromagnetic theories. 

4.1 Conjugate Spectral Intensities 
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The radiative transfer equation Eq.{l) does not remain invariant under the inverse 
transformation of the propagation vector J. Thus, it is appropriate to define two indepen- 
dent spectral intensities; /+ and I~ which satisfy the following pair of equations. 

1 dl^ 

(25) 

The equation governing is the radiative transfer equation with the propagation 

vector being equal to 5, whereas the conjugate equation governing I~ has a propagation 

vector of ~~s. 

Following a mathematical proof similar to that which was presented in section 2.1, 
one can show that the conjugate spectral intensity vectors are irrotational and therefore 
derivable from two conjugate radiation potentials i.e. 

(26) 

Substitution of Eq.(26) into (25) gives 

V • (27) 

’'■here are first order vectorial operators 

r* = ~±V + (3,j (28) 

Two vectors can be decomposed into two components as follows 

= Vi/;* -f. V X ly* (29) 

Where and W * are scalar and vector potentials respectively. Two conjugate scalar 
potentials i/>* are governed by the following Poisson equation 

V^0^ = (30) 
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whereas the vectorial potential functions obey the following homogeneous vectorial equa 
tion with repeated curl operators. 


V X V X = 0 ( 32 ^ 

The absence of the source term in Eq.(31) is demonstrated, initially, by applying a 
curl operator in Eq.(29). The result is 

V X = VxVx»-* = Vx(f|±v^ ,32, 

Secondly, the irrotationality of is proved by showing that the first and third terms 

appearing on the right hand side of Eq.(32) vanish. The proof is shown in the following 

V X = V X (J I /±3' . df-) = , X /i, = 0 (33) 

In order to establish the structural equivalency between the radiative transfer and 
electromagnetic theory, one has to demonstrate the similarity in the type of the partial 
differential operations for the conjugate radiation potentials and the electromagnetic p<v 
lentials. This is discussed in the next subsection. 


4.2 Conjugate Radiative Wave Equations: Klein-Gorden Equations 
According to the irrotational characteristcis of one expresses ' 

By applying the vectorial conjugate operator to Eq.(34), one obtains 

= L^ 


12 



The above equations, (35), are rewritten into two components explicitly as follows 


1 d^4>t 

c2 


- + 2 


0r 94>t 

c dt 


1 d^4>Z 
C2 at 2 


- vV; +2 


d<i>u 

c dt 


+ {Pi - 3 ■ ^0u)<i>t = ( - ^ + 0.S) ■ V0:r - PuS, (36) 

c at 

+ {Pl-\-S-Vp^)<i,: =[-^+P^S)-V^i^: +PuS, (37) 

c at 


where are given by 






5n f - f. 


(38) 


in which 0*(^o) potentials and their normal derivatives on the 

boundary surface respectively. 


The pair of the wave equations (36) and (37) are Klein-Gordon Equations with a 
friction force represented by the term involving the first order derivative with respect to 
time. The difference between two conjugate wave equations (36) and (37) is the sign of the 
second restoring term Pi ^ 3 ■ VPi, and the source terms appear on the, ±PvS„, present 
on the right hand side of Eqs.(36) and (37). The structural similarity between the two 
theories; RT and EM is the key feature of the present theory. Presently, however, the 
functional interrelation between the properties of the two fields has not being indenti- 
fied. Thus the complete similarity between two fields can not be fully established. It may 
be mentioned, nevertheless, that there is an indirect method for the interrelation of the 
energy-stress tensors of an electromagnetic field and that of a photon gas field. Such photon 
gas-electromagnetic equivalence® has been established in the frame work of the relativis- 
tic ally covariant theory of photon gas conservation laws derived from the conservation of 
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proper photon density. This latter equation corresponds to the covariant radiative transfer 
equation in a four-dimensional space. 

Observation of Eqs.(36) and (37) reveals the following basic features. First, two 
equations are reduced to a conventional wave equation when vanishes. This is consistent 
With the well known result of the propogation of an electromagnetic wave in a vacuum. 
Secondly, in an optically thin medium, i.e., 0, « i restoring forces for two conjugate 
waves are different from each other depending on the sign of the gradient of 0^ in the 
direction of the propagation. For example, the coefficient of a restoring force is positive 
in the direction of the increasing 0^. However in an optically thick media /?,,>>!, the 
restoring force coefficient is positive for two conjugate waves. Thirdly, an opacity retards 
the wave propagation. Finally, the radiation potential is degenerated by the sink 0^S^ 
whereas it is generated by a source term 0,3^ for the wave. Furthermore both 
waves are partially generated by the source term proportional to the temporal variation of 

the gradient of an associated potential and the product of the opacity with the latter 
quantity. 

5. CONCLUSION 

The potential theory of radiation presented in this paper provides a new physical 
perception of the structure of a radiation field through the identification of the basic 
characteristics of the irrotationality of a spectral intensity vector. Scalar potential repre- 
.sentation of a radiation field allows a unique interpretation of the laws of radiative decay, 
effects of the boundary values on the field variables for the case of steady state, propa- 
gation, retardation, as well as the generation of a radiation wave for unsteady radiation. 
This preliminary theoretical development is focused ,by design, upon a narrowly scoped 
mathematical problem, and on the identification of the basic commonality between radia- 


live transfer and electromagnetic theories. The comparative study yields a confirmation of 
a structural similarity between two fields; radiation and electromagnetic . It is projected, 
without proof, that the ultimate equivalence theory will serve to determine the validity 
limit of the existing phenomenogical theory of radiative transfer. The results of such a 
comparative study will serve to guide the formulation of a new radiative transfer theory. 

In the mean time, the steady state potential equation seems to be a viable equation 
for the prediction of radiation coupled problems because of its accuracy and elliptic type 
ot equation. This formulation has not presently been applied in modern CFD calculations. 
However, in view of an increasing interest in radiation phenomena in high performance 
liquid rocket engines and hypersonic flows, the present method is recommended for further 
advancement toward the numerical prediction of the radiation effects on a gas phase flow as 
well as for spray combustion processes. This potential method is also expected to provide 
a unique mathematical tool for the prediction of the radiation fields of various geometrical 
objects. Such options should be explored to enhance the radiation modelling activity. 

Finally, the unsteady radiation phenomena is a largely unexplored area. The effects 
ot highly pulsed radiation on propellent heating, vaporization and ignition have a practical 
significance on the engine stajt-up and instability of liquid rocket engines. 
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